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Abstract: 

We propose entanglement for hybrid optomechanical system consisting of Bose-Einstein conden- 
sate (BEC) inside a single-mode high-Q Fabry-Perot cavity with a vibrating end mirror (mechanical 
mirror). The intracavity field couples the vibrating end mirror with collective atomic density of the 
BEC. We show that the radiation pressure generates the stationary entanglement of three bipartite 
subsystems, i.e, field-mechanical mirror, field-BEC and mechanical mirror-atoms. The resulting 
entanglement is fragile with respect to temperature. 
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Nano-optomechanical systems coupled with Bose- 
Einstein condensate provide an interesting play-ground 
both in theory and experiment. These systems are 
promising in quantum informatics [1] and quantum 
metrology [2]. Quantum correlations in these systems 
[3] generates multipartite entanglement [4], thus these 
arc suitable to study quantum telcportation [5], quan- 
tum telecloning [6] and entanglement swapping [7]. The 
most challenging goals of the modern experimental quan- 
tum mechanics is to create multi-partite entangled states 
[4]. Quantum entanglement has been widely studied in 
different systems such as optomechanical systems [8, 9], 
two atomic ensembles by sending pulses of coherent light 
through two atomic vapor cells [10], Bosc-Einstcin con- 
densates trapped in double well [11] and in optical lattice 
[12]. In this paper we develop steady state off-resonant 
multipartite entanglement in hybrid system formed by 
Bose-Einstein condensate (BEC) inside a cavity with 
high-finesse, single mode optical cavity, and a moving 
end-mirror. The atom-light interaction is enhanced as 
the condensate atoms are collectively coupled to the same 
light mode. The intracavity field acts as a nonlinear 
spring which couples the BEC atoms with vibrating mir- 
ror of the cavity. We measure the entanglement between 
mechanical and atomic modes with intracavity field and 
also between mechanical and atomic mode themselves. 
The entanglement capabilities of optomechanical system 
arc modified due to the back action induced by the 
atoms. We show that the entanglement is sensitive to 
temperature variation. Though bipartite entanglement, 
namely, atom-field and mirror-field entanglement survive 
at higher temperatures, however atom-mirror entangle- 
ment is relatively fragile and available at lower tempera- 
ture scale. The use of mutual coupling between mechan- 
ical and atomic subsystems provides coherent quantum 
control at mesoscopic scale [13]. 

We consider a Bose-Einstein condensate (BEC) of N 
two-level atoms strongly interacting with a quantized sin- 
gle cavity mode of frequency, uj c . The field inside the cav- 
ity forms one-dimensional optical lattice potential. The 
intracavity field is coupled to the external fields incident 



from one side mirror with partial reflectivity, when the 
cavity is coherently driven by a laser light with frequency, 
wl, and amplitude, E. We consider the other end mirror, 
with 100% reflectivity, of the optical cavity of length, 
L, as moving and following harmonic oscillation. The 
mirror oscillates with frequency, u) m , and in the absence 
of the radiation-pressure coupling it undergoes Brown- 
ian motion as it is connected with thermal environment. 
The system is open as cavity field is damped due to the 
leakage of photons through the fixed mirror, and the vi- 
brating end mirror is connected to a bath, at finite tem- 
perature T. The Hamiltonian of the system made out of 
the intracavity field, the BEC, and the moving end mirror 
of the cavity is 



H = H n 



(1) 



where, H m describes moving-end-mirror and its coupling 
to the light field, H a describes BEC and its coupling to 
the intracavity field and ifr accounts for dissipation and 
coupling of subsystems to the thermal reservoirs. 

The mirror-field Hamiltonian H m is given explicitly 
[14] as, 
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(2) 



where, A c = w c - wl- Dimensionlcss momentum and 
position operators of the mechanical oscillator of mass 
m are, respectively, p m and q m , with commutation re- 
lation [q m ,Pm] = i- Furthermore, c and c' are the 
annihilation and creation operators, for the intracavity 
field which satisfy the commutation relation [c, c^] = 1. 
The input laser field of amplitude |E| = \J P 'n / 'huj-^, 
populates the cavity mode which is coupled to the me- 
chanical oscillator (vibrating mirror) with frequency w m 
through radiation-pressure via the coupling parameter, 
Cmc = u> c \/fi/muj m /L. 

In order to describe the motion of the BEC atoms in- 
side an optomechanical cavity, we assume that the atoms 
are trapped in one-dimensional optical lattice and their 
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motion is quantized along the cavity axis. We also as- 
sume that the atom-field detuning A a = cjl — w a is large, 
so that spontaneous emission is negligible, and we can 
adiabatically eliminate the internal excited state dynam- 
ics of the atoms. Our analysis is valid for weakly inter- 
acting BECs with little or no interactions, a situation 
that can be realized experimentally with Feshbach scat- 
tering resonances [15]. In addition, following the discus- 
sion in Ref. [16, 17], we consider that the analysis is also 
valid for strongly interacting homogeneous condensate, 
where, nonlinear term can be replaced by an effective 
potential provided the external modulation causes slight 
changes in density profile of the condensate. In the pres- 
ence of atom atom interaction, however, we employ Bosc- 
Hubbard Hamiltonian [18]. The atom- field Hamiltonian, 
H a , is written as [19], 



2 m a dx 2 



n — t — c'c 



$(x)dx, (3) 



where, is the bosonic field annihilation operator 

for the atoms. The intracavity mode couples to the 
BEC atoms through the dipole interaction via g{x) = 
g a cos(kx), where, k is the wave number of the light field. 
The cavity dynamics of BEC can be described in a homo- 
geneous two-mode model where, the macroscopically zero 
momentum state is only coupled to symmetric momen- 
tum states ±2fifc via absorption and stimulated emission 
of cavity photons [20]. Accordingly, we can write the 
atomic field operator ^(x) as, 



ii(x) = [b + V2cos(2kx) S 2 ]/vT, 



(4) 



where, b and 62 being the bosonic annihilation operator 
of the corresponding modes. By inserting the ansatz into 
Eq.(3), we get the following second quantized Hamilto- 
nian operator in Bogoliubov approximation 



HU N , fi^,„ 2 „ 2 N , 

H a = - cJc+ — (p a + q a ) + Cacfic l cg a , 



(5) 



where, il = 4w r = 2hk 2 /m a . The parameter U = gil A a 
is the optical lattice barrier depth per photon and rep- 
resents the atomic back action on the field [21]. Here 
Cac = V / NC/ /2, i.e., the single-atom interaction U Q is 
increased by the square root of number of the atoms, 
N. Along the cavity axis (x-axis), the intracavity field 
forms an optical lattice potential of period A/2, and 
depth, HU &c. The Eq.(5) describes a mechanical os- 
cillator coupled to the cavity field through the radia- 
tion pressure, where position and momentum operators, 
q a = (b + b')/\^2 and p a = (b — b')/iy/2, satisfy the com- 
mutation relation [q a ,p a ] = i. We call this fictitious mir- 
ror as atomic-mirror, which is an analogy to the moving- 
end-mirror of the cavity. 

For the full description of the hybrid optomecahni- 
cal system, we must include the effects of the dissipation, 
on the intracavity field, the damping of the mechanical 
oscillator, and the damping of the atomic-mirror. The 



Hamiltonian Ht in Eq.(l) accounts for these processes, 
and these noise processes are included via standard quan- 
tum noise operators [22]. 

In order to describe the complete dynamics of the 
subsystems involved, an adequate choice is to use the 
formalism of the quantum Langcvin equations. The ex- 
plicit form of the set of Langevin equations for the general 
system at hand reads as, 



c = (-iA + i( mc q m - i( ac q a - k)c + E - 

pm = -U> m q m + (mc C f C - 7 m p m + / m , 

q a = ttp a , 

p a = -£lq a - Cac c f c, 



Qihi Cii 



(0) 



where, A Q = A c + NU /2. For simplicity we omit the 
hat sign from the operators in Eq.(6) and in later cal- 
culations. Here, n and 7 m respectively characterize the 
dissipation of the cavity field, mechanical oscillator and 
collectively density excitations of the BEC, respectively. 
The cavity input noise is delta correlated in time do- 
main, i.e, (c in (t) c} a (t')) = S(t -t') for Hlu c /K b T » 1 
and all other correlations are zero. Mechanical Brownian 
noise operator / m with zero mean value is generally non- 
Markovian. However, the mechanical frequency never 
becomes larger then hundreds of MHz and even for cryo- 
genic temperature, the correlation function of / m can be 
approximated as (f nl {t)f m (t r )) = 7 m (2n+l)S(t-t') [23], 
where, n = [cxp {fux) m / K^T} — 1] _1 is the equilibrium 
phonon number of the mechanical oscillator. 

We now rewrite each Hciscnberg operator of the 
Eq.(6) as the sum of its steady state mean value and 
an fluctuation operator with zero mean value i.e q m = 

Qms + Sqm, Pm = Pms + $Pm, 9a = <7as + Sq a , p a = p as + Sp a , 

c = c s + Sc. Neglecting the atomic losses due to heating 
we obtain the following linearized Hcisenbcrg-Langevin 
equations 



Sq m = u m Spm , 

Spm = -W m Sq m + Xmc SX - 7 m Sp m + f m , 

Sq a = ilSp a , 

5p a = -VL5q a - XacdX , 

SX = A SY — k SX + \/2~kX VCl , 

SY = -ASX + XmcSq m - XncSQ - kSY 



'2kY w 



where, A = A Q - |c s | 2 (C mc /u; m - Q c /Si) is the effective 
cavity detuning . Here we consider the phase of the cav- 
ity field so that c s is real and positive, and we have de- 
fined the cavity field quadratures SX = (Sc + Sc { )/^/2 
and SY = (Sc — 5c^)/iy/2, and the corresponding Her- 
mitian input noise operators 5X m = (Sc- m + 5c\ n )/\/2 
and SY[ n = (Sc — Scl n )/iy/2. The linearized quantum 
Langevin equations show that the fluctuations of me- 
chanical mirror and atomic-mirror are now coupled to 
the cavity field quadrature fluctuations by the effective 
couplings Xmc = Cmc c s %/2 and Xac = Cac c s V2, which can 
be made very large by increasing the amplitude c s of the 
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intracavity field. Moreover, significant entanglement of 
optical field with mechanical oscillator and atomic mir- 
ror become possible. 

The system of linearized quantum Langevin Eqs.(7) 
can be written in compact matrix form as R(t) = 
M R(t) + F(t), where, R is the vector of the quadrature 
fluctuations and F correspond to noises. Here M is drift 
matrix. Since the quantum noises are white in nature and 
the dynamics is linearized, the state of the system will be 
continuous variable (CV) tripartite Gaussian state, and 
is completely determined by first and second moments. 
The equations of motions for first and second moments 
are 



(8) 



d = Md, 

V = MV + VM t + V F , 

where, d = (R) is the displacement vector and Vij = 
((Ri Rj + Rj Ri) — 2di dj) is 6 x 6 covariance matrix. Here 
Vp is the covariance matrix of noises and in superscript t 
describes transpose of matrix. The system will be stable 
and reaches its steady state only if all the eigenvalues of 
the drift matrix M are in the left half plane. The stability 
conditions can be obtained by applying Routh-Hurwitz 
criterion [24]. In steady state covariance matrix fulfills 
the Lyapunove equation, 



MV + VM l 



-Vp 



(9) 



which is the linear matrix equation and can be straight 
forwardly solved however the general exact expression is 
too cumbersome and is not reported here. One can ex- 
tract all the information about the steady state of the 
system from the correlation matrix. The entanglement 
between any two of the three bipartite states can be mea- 
sured by tracing out one of the three modes, that is, me- 
chanical, atomic or optical mode. 

We can quantify the steady state entanglement by 
considering the logarithmic negativity, Em [25]. As Em 
only measures the entanglement of bipartite system, we 
measure the entanglement of atomic-mirror and moving 
end mirror separately with optical field, in addition we 
also measure the entanglement between atomic-mirror 
and moving end mirror. We denote the logarithmic neg- 
ativities for the mechanical mirror-field, atomic mirror- 
field and mechanical mirror atomic mirror bimodal par- 
tition as E mc , E ac and £ ma , respectively. 

At first we measure the entanglement between me- 
chanical mirror and optical field E mc , which is obtained 
by tracing out the atomic- mirror mode, i.e, removing the 
rows and columns of V correspond to atomic-mirror. The 
reduced state is still Gaussian and fully characterized by 
4x4 matrix V mc and is given by 



V m 



X z 
z l y 



Here, X, y, Z are 2x2 matrices, where matrix X corre- 
sponds to the mechanical mirror, y refer to the intracav- 
ity mode and Z describes the correlation between these 



two systems. In order to measure the entanglement be- 
tween mechanical mirror and optical field, we consider 
the logarithmic negativity E mc , which, in case of CV, 
E mc can be defined as [26] 



E mc = max[0, — In 2 e] , 



(10) 



where, e = 2- 1 /2 E ( Klc )-E(K lc ) 2 -4detK lc ] 1 /2}i/2, 
with XX^mc) = det X + dety - 2detZ. The Gaussian 
state gets entangled only if the eigenvalue e < 1/2, and 
it is Simon's necessary and sufficient entanglement non- 
positive partial transpose criterion of the Gaussian states 
[27]. This condition can also be written as 4detV mc < 

E(^mc) - 1/4. 

For our numerical calculation we take the parame- 
ters from the experimental work reported in [13, 20, 28]. 
In Fig.l we plot the steady state entanglement between 
mechanical mirror and optical field Fig. 1(a), atomic mir- 
ror and optical field Fig. 1(b) and mechanical and atomic 
mirrors Fig. 1(c) measured by the logarithmic negativ- 
ities -Emc, E^c and E ma , respectively. In Fig. 1(a) we 
plot the E mc against normalized detuning A/uj m and 
coupling rate of the mechanical mirror with cavity field 
£ mc for £ ac = 0.7 e mc and co m ~ Q. The entanglement 
measure between mechanical mirror and optical field in- 
creases with £ mc and is present only for small interval of 
values of A around A ~ 0.5 uj m . 

Similarly, the steady state entanglement of the 
atomic mirror and intracavity field is measured by the 
logarithmic negativity E ac as a function of the normal- 
ized detuning and the coupling of the atomic mode with 
intracavity field for ( = 0.01 £ mc and oj m ~ f2, as shown 
in Fig. 1(b). The features observed in the entanglement 
measure between atomic and cavity modes are similar to 
those present in E mc . The Entanglement measure be- 
tween mechanical mirror and optical field increases with 
£ ac and present only for small interval of values of A 
around A ~ 0.5 w m . 

Analogously, we measure the entanglement between 
mechanical mirror and atomic mirror. We numerically 
calculate the logarithmic negativity, E^, associated with 
the steady state correlation matrix formed by the atomic 
and mechanical modes, by tracing out the cavity mode 
shown in Fig. 1(c). Here the logarithmic negativity, E ma , 
is calculated as a function of the normalized coupling of 
the mechanical ( mc /ui m and atomic mirrors with intra- 
cavity field Cac/^m for A ~ 0.6 w m and w m ~ fl. 

It is important to understand the behavior of entan- 
glement with respect to the temperature T. In Fig. 2 we 
show the logarithmic negativity Ei, where, (i = mc, ac, 
ma) as a function of the environmental temperature T. 
In Fig. 2(a), it is noted that the entanglement between 
mechanical mirror and intracavity field is very robust 
with respect to temperature. The black solid line for 
Cac = shows a significant amount of entanglement is 
present upto temperature 60 K. The red dashed line for 
Cac = 0.4 C mc in Fig. 2(a) shows that the entanglement 
between mechanical mirror and intracavity mode persist 
only upto 20 K. Similarly, the robustness of the entan- 
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(a) (b) (c) 

FIG. 1: (Color online) (a) Plot of logarithmic negativity E mc as a function of the normalized detuning A/u m and coupling of 
mechanical mirror with optical field £ mc for Cac = 0.7 Cmc and u m ~ fi. The optical cavity length L = 1mm, driven by a laser 
with wavelength A = 1000 nm, power P = 50 mW and u) m — fi. The mechanical mirror has a frequency LOm/2n = 10 MHz, 
and damping rate ■y rn /2ir = 100 Hz, its temperature is T = 100 mK and cavity finesse is F — 1.07 x 10 4 .(b) _B ac against A/ujm 
and Cac for Cmc = i£ = iuj c yh/mUm./^j (i = 0.01) and U) m — Q,. (c) _E ma against normalized coupling of mechanical mirror 
Cmc = Cmc/t^m and atomic mirror Cac = Cac/ W m for A = 0.6 cj m and temperature is T = 1 /iK. 




FIG. 2: (Color online) (a) Plot of logarithmic negativity _E mc as a function of environmental temperature T for A ~ 0.5 cj m 
and com — fi. The black line (solid) for Cac = and the red line (dashed) for Cac = 0.4 Cmc -(b) _E ac against the environmental 
temperature T for A ~ 0.5 U) m , Cac = 100 Hz and ui m ~ Q. Cmc = if = iu} c ^/h/muj m /L, the black (solid) line for i = and red 
dashed line (dashed) for i = 0.4 (c) -E ma against temperature T for A ~ 0.6 aj m and ui m — fi. The black (solid) and red dashed 
(dashed) lines for Cmc = Cac = aj m and Cmc = Cac = 0.001 aj m respectively. The other parameters are the same as in Fig.l. 



glement between atomic mirror and intracavity mode as 
a function of the reservoir temperature of the mechan- 
ical mirror is shown in Fig. 2(b). The black solid and 
red dashed lines are for Cmc — 0.1 C and Cmc = 0.4 C, re- 
spectively. The BEC mode absorbs some phonons taken 
from the mechanical mirror by the field and its temper- 
ature is increased. There is no direct connection be- 
tween BEC and thermal environment for zero interac- 
tion between mechanical mirror and cavity field, hence 
E ac is independent of temperature as the interaction is 
zero, i.e, Cmc = 0. In Fig. 2(b), the black solid line 
for a smaller value Cmc = 0.1 C, of the interaction be- 
tween mechanical mirror and intracavity field, shows that 
E ac persists upto 20 K. The red dashed line in Fig. 2(b) 
shows the entanglement between atomic and cavity field 
modes only presents upto 5K as the coupling between 
mechanical mirror and intracavity is increased, that is, 



Cmc = 0.4 C- The steady state entanglement between me- 
chanical mode and atomic mode as a function of the tem- 
perature is shown in Fig. 2(c) for two different values of 
the coupling of the cavity field with atomic and mechan- 
ical modes. The E ma is fragile with respect to temper- 
ature. The black solid line and red dashed line are for 
Cmc = Cac = oj m and Cmc = Cac = 0.001 u m , respectively. 

In Fig. 3 we plot the steady state entanglement of me- 
chanical mirror with intracavity field (black line) versus 
normalized detuning A/w m and compere it with, steady 
state entanglement of the BEC mode with intracavity 
mode (red dashed line) and the steady state entangle- 
ment of the BEC mode and mechanical mode (green dot- 
ted line). Fig. 3 shows that the mirror-field and atom- field 
entanglement exist simultaneously, whereas atom-mirror 
entanglement is observed only for large values of the de- 
tuning. 
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FIG. 3: (color online) Plot of Logarithmic negativity En as 
function of normalized detuning A/cj m for cj m = 10 MHz and 
SI = 0.9 w m . Mirror- field entanglement in the absence of BEC 
mode (black solid line) and BEC-field entanglement in the 
absence of mechanical mode (red dashed line). Moreover, the 
entanglement between mechanical mode and BEC mode in 
the absence of cavity mode (green dotted line) for = u m — 
2tt x 1MHz. Furthermore, ( ac = 200 Hz and Cmc = 300 Hz 
and temperature T = 1 /iK. The other parameters are the 
same as in Fig.l. 



For experimental realization of the generated entan- 
glement, one has to measure several quadrature correla- 
tions [29] and these have been experimentally measured 
for the entanglement of the two optical modes [30] . How- 
ever, in our case we consider another Fabry Perot cavity 
adjacent to first one and is driven by a weak laser field, 
a scheme of this kind has been discussed in [9], in which 
the position and momentum of the mechanical mirror is 
directly measured by homodyning the out put of the sec- 
ond cavity for suitable values of the detuning and decay 
rate of the second cavity whereas the quadratures are di- 
rectly measured by homodyning out put of the cavity. 
Finally, the quadratures of the BEC mode are measured 
by homomodning a weak field fed into the cavity as in 
Ref. [31]. 

In conclusion, we have studied the hybrid quantum 
correlation of tripartite system which consists of opti- 
cal cavity, a Bose-Einstein condensate mode, and a me- 
chanical mirror mode of the Fabry-Perot cavity. We 
have shown steady state entanglement of the three bipar- 
tite subsystems: the mechanical mirror-field, the atomic 
mirror-field, and the mechanical mirror-atomic mirror by 
quantifying in term of logarithmic negativity. The result- 
ing entanglement is fragile with respect to temperature. 
The realization of such systems opens new horizons for 
the realization of the quantum memories [32] and inter- 
faces for the quantum information processing. 
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